Diffusion-Limited Aggregation on Curved Surfaces 
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We develop a general theory ol transport-limited aggregation phenomena occurring on curved sur- 
faces, based on stochastic iterated conformal maps and conformal projections to the complex plane. 
To illustrate the theory, we use stereographic projections to simulate diffusion-hmited-aggregation 
(DLA) on surfaces of constant Gaussian curvature, including the sphere [K > 0) and pseudo-sphere 
{K < 0), which approximate "bumps" and "saddles" in smooth surfaces, respectively. Although 
curvature affects the global morphology of the aggregates, the fractal dimension (in the curved 
metric) is remarkably insensitive to curvature, as long as the particle size is much smaller than the 
radius of curvature. We conjecture that all aggregates grown by conformally invariant transport 
on curved surfaces have the same fractal dimension as DLA in the plane. Our simulations suggest, 
however, that the multifractal dimensions increase from hyperbolic {K < 0) to elliptic {K > 0) 
geometry, which we attribute to curvature-dependent screening of tip branching. 

PACS numbers; 61.43.Hv, 47.54.-|-r, 89.75.Kd 



The Laplacian growth model and its stochastic ana- 
logue, diffusion-limited aggregation (DLA) [3], describe 
the essential physics of non-equihbrium pattern forma- 
tion in diverse situations [2], such as viscous finger- 
ing [7] , (quasi-static) dendrite sohdification [5] , dielectric 
breakdown [6], and dissolution [12], depending on con- 
ditions at the moving, free boundary. Some extensions 
to non-Laplacian growth phenomena, such as advection- 
diffusion-limited aggregation [14, 17] (ADLA) and brit- 
tle fracture [13], are also available. Almost all prior work 
with these models has assumed flat Euclidean space, typ- 
icahy a two-dimensional plane, but real aggregates, such 
as mineral dendrites [18], cell colonies [19], and cancer- 
ous tumors [20], often grow on curved or rough surfaces. 
Aside from a few studies of Eden-like clusters [19] and 
continuous viscous fingers [21, 22] on spheres, it seems 
that the effects of surface curvature on pattern forma- 
tion have not been investigated. 

In this Letter, we extend transport-limited growth 
models to curved two-dimensional surfaces via confor- 
mal projections from the plane. Time-dependent confor- 
mal maps are widely used in physics [1] and materials 
science [15] to describe interfacial dynamics in two di- 
mensions. Continuous conformal maps have long been 
applied to viscous fingering [7, 8], and more recently, 
Hastings and Levitov introduced stochastic, iterated con- 
formal maps for DLA [25]. Both continuous and stochas- 
tic conformal-map dynamics have also been extended to 
other conformally invariant (but non-Laplacian and non- 
linear) gradient-driven transport processes [10, 14, 17], 
such as advection-diffusion in a potential flow [9, 11] or 
electrochemical transport in a quasi-neutral solution [10]. 
Indeed, there is nothing special about harmonic functions 
(solutions to Laplace's equation) in the plane, aside from 
the direct connection to analytic functions of a complex 
variable (real or imaginary part). The key property of 
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FIG. 1: Stereographic projection, $~ , from the exterior of 
the growing object f2m(t) on a sphere of radius R to the exte- 
rior of the shadow, Qz{t), on a complex plane. The point $(2) 
is projected from the north pole to the point z. The origin 
of the 2;-plane is tangent to the sphere at the south pole, and 
the latitudinal angle </) is measured from the south pole. 



conformal invariance is shared by other equations [10, 15] 
and, as note here, applies equally well to conformal 
(i.e. angle preserving) transformations between curved 
surfaces. Here, we formulate continuous and discrete 
conformal-map dynamics for surfaces of constant curva- 
ture, by a sequence of mappings from the complex plane, 
and we use the approach to study the fractal and multi- 
fractal properties of DLA on curved surfaces. 

Transport-limited growth on curved surfaces. — At 
first, it would seem that conformal-map dynamics can- 
not be directly applied to a non-Euclidean geometry, 
and certainly the formulation could not be based only 
upon analytic functions of the complex plane. Nev- 
ertheless, if there exists a conformal (angle-preserving) 
map between the curved manifold and the plane, and 
if the underlying transport process is conformally invari- 
ant [10, 11, 14, 15, 17], then any solution to the transport 
equations in the non-Euclidean geometry can be confor- 
mally mapped to a solution in the complex plane, without 



2 



(a) 



(b) 





South Pole 

FIG. 2: DLA clusters on the elliptic (left) and the hyper- 
bolic (right) geometries. The elliptic geometry is isometri- 
cally embedded on the surface of a sphere and the hyperbolic 
geometry is visualized on the Poincare disk with the metric, 
ds = dz/{l- {\z\/2Rf). We use R/y/To = 100 and aggre- 
gate 4942 and 9001 particles to fill the great circles of radius 
{n/2)R (dashed lines) in the elliptic and hyperbolic geome- 
tries respectively. 



changing its functional form. Such maps do exist, as we 
illustrate below with important special cases. 

Let ^rn{t) be the exterior of a growing object on a 
non-Euclidean manifold, M, and $ be a conformal map 
from a (part of) complex plane to M. We speak of the 
domain ^lm{t) as having its shadow, ftzit) = ^~^{Qm(t)), 
on the complex plane under the projection . As in the 
flat surface case, we can describe the growth by a time- 
dependent conformal map, g{w,t), from the exterior of 
the unit disk, D,^, to the exterior of the growing shadow, 
i^z{t)^ which in turn is mapped onto the curved geometry 
by the inverse projection. Care must only be taken that 
the dynamics of g{w,t) should describe flz{t) in such a 
way that the evolving object, flra{t), follows the correct 
physics of growth on M, rather than in the intermediate 
complex plane, which is purely a mathematical construct. 

Let us illustrate this point for both continuous and dis- 
crete versions of conformally invariant, transport-limited 
growth [15, 17]. For continuous growth, we generalize 
the Polubarinova-Galin equation equation [15-17] for a 
curved manifold with a conformal map ^{z) as follows 



Re{w g'{w) gt{w)} 



a a{w, t) 
Wog{w)\^ 



(1) 



for \w\ = 1, where a is a constant and (t(w, t) is the 
time-dependent flux density on the boundary of VL^ ■ This 
result is easily obtained by substituting ^ o g for g in the 
original equation. 

For stochastic growth, we adjust the Hastings-Levitov 
algorithm [25] on the shadow domain. The algorithm is 
based on the recursive updates of the map, 

gn[w) ^ gn-i° (t)\„,e„{w), gn{w) ^ g{w,tn), (2) 

where (j)x^g is a specific map that slightly distort fl^ by 
a bump of area A around the angle 9. While the ran- 
dom sequence {On} follows the probability distribution, 
p{9, tn) (x cr{e^^ , tn), invariant under conformal maps, the 



preimage of bump area, A„, should be determined so that 
the bump area is fixed as Aq on the manifold, M. Thus 
the bump size of the n-th accretion is determined by 
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(3) 



We are not aware of any prior modification of the 
Hastings-Levitov algorithm with this general form. Pre- 
vious studies on DLA in a channel geometry [28] can be 
viewed as an example with M — {z : < arg z < 27r} 
and $(z) = log(2;) although the manifold is Euclidean. 

Stereographic projections. To illustrate the general 
theory, we first make use of classical stereographic pro- 
jection [27] to describe growth on a sphere. Stereo- 
graphic projection is obtained by projecting the surface 
of a sphere from the north pole to a plane whose origin 
is tangent to the south pole; see Fig. 1. If $ is an in- 
verse stereographic projection with sphere of radius R, 
maps the point {R, 0, 9) in spherical coordinates to 
z = i?tan(0/2)e'^ in the complex plane. Here is the 
azimuthal angle and is the latitudinal angle measured 
from the south pole. If the modulus, ] • ], on the sphere 
is defined to be distance to the origin (south pole) in the 
curved metric (arc-length of the great circle) in a similar 
way to ] • ) on a complex plane, \z\ and 1$(2:)1 satisfy 
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(4) 



The Jacobian factor of the projection is angle- 
independent; thus, it is given by 



l$'(z)l = 



dmz)\ 

d\z\ 



1 



l + i\z\/2Ry 



(5) 



from the derivative of Eq. (4). Now Eq. (5) can be used to 
continuous dynamics, Eq. (1), and stochastic dynamics, 
Eq. (3) on sphere. 

While the surface of a sphere is a three-dimensional vi- 
sualization of elliptic (or Riemannian) geometry, we can 
also obtain a conformal projection from hyperbolic geom- 
etry to a complex plane as well. Unlike elliptic geometry, 
hyperbolic geometry can not be isometrically embedded 
into 3D Euclidean space; only a part of the geometry 
can be embedded into 3D as a surface known as pseu- 
dosphere. Hyperbolic geometry has a negative constant 
curvature, K = — 1/i?^, as opposed to the positive one, 
K = l/R^, of elliptic geometry. The projection can be 
obtained by simply viewing hyperbolic geometry as a sur- 
face of a sphere with an imaginary radius, iR, as sug- 
gested from the sign of curvature. The projection can be 
still defined; substituting iR for R alters Eqs. (4) and (5) 
as 



2R \ 2R 



and 



l$'(z) 



1 



\z\/2Rr- 
(6) 
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Geometry 


Elliptic 


Euclidean 


Hyperbolic 


Df 


1.704±0.001 


1.704±0.001 


1.693±0.001 


D2/D1 


0.577±0.003 


0.529±0.006 


0.527±0.001 


Ds/Df 


0.591±0.006 


0.503±0.005 


0.499±0.001 


Di/Df 


0.617±0.012 


0.486±0.004 


0.482±0.001 


D5/Df 


0.631±0.018 


0.473±0.004 


0.469±0.001 



TABLE I: The fractal dimension Df and the multifractal di- 
mensions D2q+i of DLA clusters on three different geometries. 
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FIG. 3; (a) The average number of particles, (A'^) , to fill the 
disk of conformal radius, Ai — {tt/2)R, versus radius R on 
the three geometries, (b) Moments (A^) of the pre-image 
bump size, Ajv, versus R, which define O213+1 via (9. 



The image of hyperbolic geometry under the projection 
is thus limited to the inside of a disk with radius 2R, 
and the boundary, \z\ — 2R, corresponds to the infinity. 
The stereographic projection serves as a non-isometric 
visualization of the geometry known as Poincare disk. 
The length element dz at dz in the disk corresponds to 
the element c?z/(l — (|z|/2i?)^) in the hyperbolic universe. 

DLA on constant- curvature surfaces. — Using the har- 
monic probability measure, p{9) ~ I/Stt, and the mod- 
ified Hastings-Levitov algorithm with Eq. (3), we grow 
DLA clusters on curved surfaces. Fig. 2 shows clusters 
on elliptic and hyperbolic geometries. On the Poincare 
disk, the size of particles becomes smaller and smaller as 
they approach the infinity, \z\ ^ 2R. 

An analytic advantage of the conformal mapping for- 
mulation is that the Laurent expansion of g{w, t) gives us 



information on the moments of the cluster; the confor- 
mal radius, Ai, and the center of charge, Aq, come from 
the first two terms of the expansion, g{'w) ~ Aiw + Aq. 
Such coefficients from DLA clusters on curved surfaces 
are not useful as they are the moments of the shadow, 
not the original cluster. Using the property that circles 
are mapped to circles under an (inverse) stereographic 
projection, it is possible to define analogous quantities 
corresponding to Ai and Aq. We note that the image, 
<^{Aiw -\- Aq), of the unit circle, \w\ = 1, is also a circle 
on M, and it is the one that best approximate the cluster 
and the far-field in Vlmit). Thus, we define the conformal 
radius, Ai, and the center of charge, Aq, on M to be the 
radius and the deviation of ^{Aiw -\- Aq) respectively: 



A = ^ { MA, + l^ol 



MA, - l^oDI }, (7) 
MA, - \A,\)\ }. (8) 



On curved surfaces, however, the fractal dimension Df, 
cannot be determined from the scaling, n ~ (^i)^^ , since 
the geometry is not linearly scalable with A, ; the two do- 
mains within radii of different values of A, are not self- 
similar to each other, and the log-log plot of n versus A, 
is not linear either. Instead, we use the sphere radius R 
as a relevant length scale while keeping Ai proportional 
to R for self-similarity. Thus, we grow a cluster until 
Ai reaches at (j)oR for various radius R, but with a fixed 
particle size, Ao — 1, and a fixed angle (t>o = 7r/2. The 
angle (/jq is as such since the circle (dashed lines in Fig. 2) 
becomes a great circle on a sphere. If N is the number 
of particles to fill the radius, the fractal dimension is de- 
termined from {N) ^ R^f . From the statistics of 1000 
clusters for each of i? in a geometrically increasing se- 
quence from 79 to 400, we obtain the fractal dimensions, 
Df w 1.70, as shown in Table I. Fig. 3(a) shows {N) 
versus R in the three geometries. The relative deviation 
of Df between geometries is surprisingly small compared 
to the deviation of surface properties caused by the cur- 
vature. For example, the area within the radius {tt /2)R 
on elliptic (or hyperbolic) geometry is about 23% smaller 
(or larger) than the corresponding area on the Euclidean 
surface; however, this factor only change the pref actor, 
not the exponent, of the scaling, (A'') ~ R^f . 

Our results suggest that Df is insensitive to the cur- 
vature because, on small length scales comparable to the 
the particle size \/Ao, the surface is locally Euclidean. 
This result is consistent with ADLA [14, 17] whose frac- 
tal dimension is not affected by the background fluid flow. 
We conjecture that any conformally invariant transport- 
limited aggregation process on a d-dimensional curved 
manifold has same, universal fractal dimension as DLA 
in flat (Euclidean) space in d dimensions. 

More subtle statistics of the aggregates, however, are 
revealed by multifractal scalings. The multifractal prop- 
erties [29, 30] related to from the probability measure. 
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however, seem to depend on the curvature. Following 
Ref. [26], we measure the multifractal dimensions Dq 
from the relation, 



for any conformally invariant transport-limited growth 
process, in any number of dimensions. 

We acknowledge helpful discussions with Ellak Somfai. 



(9) 



The averaging in Eq. (9) is made over Ajv at uniformly 
distributed angle. Fig. 3(b) shows the first three mo- 
ments of Ajv as functions of R, and Table I shows Dg/Df 
in the three geometries. We note that Dq for elliptic ge- 
ometry are yet transient, influenced by finite-size effect. 
The multifractal dimensions do not satisfy the inequal- 
ity Dq > Dq/ for q < q' [29]. Since the space itself is 
finite and the distance to infinity is bounded, the insta- 
bility in growth is pronounced in elliptic geometry. The 
normalized center-of-charge fluctuation, {\Ao\'^y^'^ /R, is 
about 0.3 for on elliptic geometry although it decreases 
as R increases. For Euclidean and hyperbolic geome- 
tries, the fluctuations are around 0.03 and 0.02 respec- 
tively. It should also be noted that the growth probability 
on sphere is not harmonic since the far-field potential, 
— log|2;|, on plane is mapped to the singular potential, 
log(7r — (f)), aroimd the north pole. 

Even with the finite size effect, we conjecture that 
multifractal dimensions increase in the order of hyper- 
bolic, Euclidean, and elliptic geometries. The justifica- 
tion is that, when the measurements are similarly made 
for smaller latitudinal angles, (po = 7r/3 and 7r/6, Dq on 
elliptic and hyperbolic geometries converge to Dq on Eu- 
clidean plane. This also support that the slight difference 
in Dq between Euclidean and hyperbolic geometries are 
not from statistical error. We believe the dependence of 
Dq on curvature is related to the depth of fjords. On el- 
liptic geometry, the circumference of a circle with radius 
4>oR is 27r sin 4>oR: shorter than 2n(j)oR on Euclidean ge- 
ometry. If we consider a pair of branches which make the 
same opening angle on the two geometries, the area be- 
tween them is more screened on elliptic than on Euclidean 
geometry. Therefore, we expect a different distribution of 
A. Due to the competition on smaller circumference, tip- 
splitting events are subdued, as observed in continuous 
growth [22], and eventually fewer branches survive. The 
opposite argument applies to hyperbolic geometry, where 
the circumference is given by 27r sinh 4>oR, inter-branch 
area is less-screened, and tip-splitting is encouraged. In 
future work, it would be interesting to analyze the spec- 
trum of A and its connection to tip splitting events. 

In summary, we have developed a mathematical the- 
ory of transport-limited growth on curved surfaces and 
applied it to DLA on two-dimensional surfaces of con- 
stant Gaussian curvature. Our simulations suggest that 
the fractal dimension of DLA clusters is universal, in- 
dependent of curvature, and depends only on the spatial 
dimension. The multifractal properties of DLA, however, 
seem to depend on curvature, since tip-splitting is related 
to the different degrees of screening the inter-branch ar- 
eas. We conjecture that these results are hold in general. 
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